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Previous experiment on Fe/Cr/Fe(001) trilayers reported a surprising observation that the in- 
terlayer exchange coupling could be enhanced drastically by the bombardment of irradiation even 
at low fluences. We propose that it is due to the resonant states in the spacer, made possible 
when the topography of both interfaces is correlated and exhibits prominent Fourier components. 
A systematic procedure is developed to handle the interface roughness and predict on how to opti- 
mize the interlayer coupling. This method can be extended to bridge the gap between theories and 
experiments in other heterojunctions. 
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I. INTRODUCTION 

The interlayer exchange coupling has been studied for 
more than twenty years[l|, 0, y, [j, [Bl, @ with applica- 
tions in phenomena such as Giant Magnetoresistance[3, 
3 (GMR) and Tunnelling Magnetoresistance (TMR). 
When it comes to estimate the effect of interface rough- 
ness (IR) on the coupling, there are few microscopic the- 
ories other than resorting to static average [3. Ifl HoL ITU . 
A recent experiment fl2j reported that the strength of the 
coupling in Fe/Cr/Fe(001) trilayers could be modified in 
a controlled manner and enhanced even at low fluences 
by the ion beam irradiation. The ion bombardment is 
known to be an excellent tool for patterning magnetic 
areas on interfaces of multilayer system without chang- 
ing the sample structure [T3|. 

This enhancement of interlayer coupling is not ex- 
pected by the conventional approach to statistically av- 
erage over the IR. Many experiments [l^, EH EH have 
been dedicated to address this surprising finding. They 
tailored the interface into exhibiting certain patterns by 
the ion beam irradiation, but so far there is no conclu- 
sive result [ill, EH- Among the many possible coupling 
mechanisms in the trilayer systems that might give rise 
to this enhancement, the authors of Ref. [l2| have esti- 
mated the percentage of having the magnetic bridges 
and concluded that their contribution to the coupling 
was minor. The spin density wave in bulk chromium has 
been known [l^, E3 to be crucial at understanding how 
the electronic structure and the topology of Fermi sur- 
face influence the magnetic properties of Fe/Cr multilay- 
ers. However, since the enhanced coupling was observed 
with a spacer of only 8 monolayers [l2|, the indirect ex- 
change is more likely to be mediated by the conventional 
Ruderman-Kittel-Kasuya-Yosida interaction^ instead of 
spin density waves which occur [13, EH only above a crit- 
ical thickness of 30 monolayers. Another possible candi- 
date for the coupling is magnetic dipole interactions [l9j. 
made possible by IR. Nevertheless, it is not likely to dom- 
inate the coupling because the analysis of STM images 
showed that the thickness fluctuations of the spacer in 
Ref.[l2| was only one monolayer. 

In this article, we present the first microscopic theory 



based on a perturbative treatment of the IR. This mod- 
ifies the wavefunction of the conduction electrons in the 
spacer and the interlayer coupling they mediate. Mainly, 
we found two additional terms to the scattering wave- 
functions due to the topography of IR besides the usual 
static average results. How and when these extra terms 
can lead to enhancement of the interlayer coupling is ad- 
dressed. 

In section II, we discuss the energy spectrum of quan- 
tum well states (QWS) in a trilayer system, which 
is shown to be broadened by the rough interface, as 
expected [20j from taking the static average. In addi- 
tion, when the two interfaces are assigned a common pe- 
riodic variation, the spectrum will get shifted upwards. 
In section III, a resonance is predicted to occur when 
the length scale of IR matches the spacer width or the 
Fermi wavelength. We propose that the latter matching 
is responsible for the enhancement of interlayer coupling 
by the ion-beam bombardment. Conclusions and discus- 
sions are arranged in section IV, in which other coupling 
mechanisms, such as the magnetic dipole interaction and 
the biquadratic coupling, and the effect of alloying at the 
interface are addressed. To preserve the conciseness of 
the main text, detailed calculations are all arranged in 
the Appendices. 



II. ROUGH INTERFACE AND CORRECTIONS 
TO QUANTUM WELL STATES 

We start by studying a simple 2-D case: an interface 
described by x = A siii(py) and assume the effective mass 
to be the same for carriers at different sides of the inter- 
face. As a plane wave comes in with momentum (k x , k y ) 
from the left, we assume the amplitude A of the rough- 
ness to be much smaller than both l/k x and 1/p. The 
fcji^l allows us to treat the roughness as a perturba- 
tion, while pA <C 1 precludes the possible emergence of 
local states at the interface, as in water waves along the 
coastline. Under this limit, the action of static averaging 
over one rough interface can be shown to be equivalent to 
the scattering of two effective smooth interfaces, the first 
and second terms in Eq.((T]). Whiles, the topography at 
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the interface creates the third and fourth terms, of which 
the strength is proportional to A and k y is shifted by p: 
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where &® n (x,y) denotes the wave function scattered by 

a smooth interface and a£ 2 are defined in Appendix A 
which contains details of the first-order perturbation cal- 
culations. 

In a trilayer system, QWS determine the interlayer ex- 
change coupling. It is crucial to study the variation of 
QWS due to the IR, which broadens the resonance peaks, 
shifts the energy spectrum, and sometimes reduces the 
number of bound states. We find that the energy shift 
becomes large when (1) 1/p is comparable to the thick- 
ness of the spacer D and (2) the topographies on these 
two interfaces are correlated, which feature could not be 
obtained by the static averaged, 0, [l(| or another ap- 
proach which combined it with the Green's function [2 1|. 
An example is shown in Fig.l, assuming both interfaces 
exhibit the same topography x = Asin(py). From Ap- 
pendix C, it is shown that the only effect of static av- 
erage is to broaden the peaks in the energy spectrum 
that represent bound states. When the first-order per- 
turbation corrections due to the IR are added, the spec- 
trum becomes shifted and the original bound states near 
the barrier edge can become unbounded. This shift has 
been observed in the photoemission experiment [20] . Note 
that both features are more pronounced at high energy 
states because the shorter wavelength of QWS renders 
them more susceptible to the topography of the inter- 
face. We shall discuss later our finding that the p as 1/D 
that causes the biggest spectrum shift often brings about 
the greatest enhancement or suppression of the interlayer 
coupling. Other details, such as the interface intermixing 
and Brillouin zone dilation [^l, seem not to be crucial for 
qualitative predictions. 

A similar feature is found in the scattering states, 
namely, interlayer coupling can be enhanced or sup- 
pressed by the IR when (1) the topographies on both 
interface are correlated and (2) p happens to be compa- 
rable to one of the wavelengths k that satisfy the condi- 
tion for the Fano resonance. Although the shift of these 
k by IR is small, they somehow collaborate to creat a 
profound effect on the coupling. 



III. INTERLAYER EXCHANGE COUPLING 

In this section, we will study the effect of IR on the 
interlayer exchange coupling for a 3-D trilayer system, 
which consists of a spacer of width D and potential V — 
sandwiched between layers L and R with potential Vl, 
Vr, respectively. Assuming the lattice structure along 
the planar directions of the interface is ideal, the coupling 
strength is a superposition of the contributions from all 




FIG. 1: Photoemission intensity is plotted as a function of the 
energy of quantum well states when both interfaces exhibit 
the same topography x = 2Asin(py). The results for smooth 
interfaces are plotted in dashed line for comparison. The gray 
line takes care of the IR by static average, while the solid 
line further considers the corrections due to the first-order 
perturbation. The parameters chosen for this plot are m* = 1, 
the depth of quantum well is Vb = 0.5 in unit of Ry/n 2 , 
2 A = 0.5 in unit of the Bohr radius, p = 0.5, broadening 
factor defined in Appendix C is set to be 1, and D = 6.8. 



the QWS in the interface Brillouin zoned |(| (IBZ). The 
existence of interfaces without the IR affects the total 
electron energies of the spacer, the shift of which equals 
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per unit area of interface. The reflective coefficients r L / R 
at the left/right interface are functions of the wave num- 
ber along the normal direction of the interface k± , which 
depends on the energy E and fcy = fc — k±. The integer n 
denotes the number of round-trip reflections by the two 
interfaces. In the presence of IR, two corrections are in- 
troduced to Eq.@: static average terms in Eq.fT]) come 
in through the r^/n, while the last two terms due to the 
perturbation contribute additional close loops when the 
topography of both interfaces are correlated. By use of 
the T-matrix and perturbation methods, wc derive the 
energy shift in Appendix B, which takes the form: 
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where d^/n are corrections after averaging the IR on 
the left/right interface, and b\/ k q {b^ q k ) are addi- 
tional transmission coefficients from momentum {k±, k y } 
to {q±,q y } (vice versa) induced by the tomography on 
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interface L (R). We checked that the above equation 
indeed reduced to Eq.([2]) at the limit of p — > 0. 

An important application of trilaycrs is the GMR and 
TMR systems where the two side layers exhibit perma- 
nent magnetization. When they are ferromagnetically 
coupled, the energy shift by Eq.Q can be calculated 
by evaluating the reflective coefficients under the corre- 
sponding Zeeman field. It is found to differ from that 
of the antiferromagnetic case. The magnetic coupling 
strength J(D) can then be defined as the difference be- 
tween these two energies @. For simplicity, we assume 
the two magnetic layers to be of the same material and 
their interface topographies are correlated. After some 
routine calculations in Appendix B, the magnetic cou- 
pling is found to be: 

n 

+ vr& (1) ir& (1) -b {1) ]j(<!±+k^D\ 

' / j I "v.k„a„ ap,k y q y ] l"p,q y k y "ap,q y k y \^ J 
q y 
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where r p /r ap denote the reflective coefficients for carriers 
with spin parallcl/antiparallel to the magnetization of the 
layer they reflect from. 

A. Uncorrelated interfaces 

An intuitive equation, which resorts to static average 
over the IR, has been widely used by experimentalists 
and theoretists alike @, i, Q3, [Tl|: 

J(D) = ^w(n)J(nd) ( 5 ) 

n 

where w(n) is the weighting of n-layers among the various 
widths bewteen two rough interfaces, d is the width of 
one unit layer in the spacer, and D is defined as the 
mean width Y] w(n)nd. This weighting w(n) is usually 
assumed to obey the Guassian distribution pi. l9l. [ToL [ill] . 

When the IR on the two interfaces are uncorrelated, 
the last term in Eq.([J| due to their interference vanishes. 
The remaining term is shown in Figf2]to be equivalent to 
the static average. It can be seen from the figure that un- 
correlated IR always reduces the coupling strength, but 
retains the same period as the exact result for smooth 
interfaces, which turns out to be roughly n/kp as pre- 
dicted by the Ruderman-Kittel-Kasuya-Yosida formula - 
a second-order result. 



B. Correlated interfaces 

When the topography of both interfaces is correlated, 
the close loops in Eq.([J]) formed by the interference will 




FIG. 2: Coupling strength in unit 10~ Ry/n eto between 
rough interfaces with random topography in 2-D is calculated 
by Eq.© (in dashed line) and the static average Eq.© (in 
gray line) as a function of the average width D in Bohr radius 
oo- The parameters we used are: the Fermi energy Ef = 1, 
V p = 0.01 and V ap = 0.1 (energy unit is Ry/i: 2 ), and the 
maximum amplitude of the IR on both interface is ao/3. For 
comparison, the result for smooth interfaces is plotted in solid 
line. 
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FIG. 3: Coupling strength is plotted as a function of the 
average film width D as Figf2]with the same parameters ex- 
cept that the two interfaces are now assumed to share the 
same topography A(y) = l/3sin(py). Again, the result for 
smooth interfaces is shown in solid line for comparison. The 
gray/dashed lines represent p = 0.05/1.5, respectively. 



survive and give rise to new signatures. The correla- 
tion can be introduced [13] by the ion irradiation through 
the trilayer, which has the capacity of rearranging the 
atoms on both interfaces. Although Ref.[l2| used an ion- 
ized helium of 5 keV, smaller than the 2-MeV Au +2 in 
Ref . [22| . we believe this correlation is still possible be- 
cause the penetration depth is inversely proportional to 
the cross section of the incoming ions. Consider the ex- 
treme case that both interfaces share the same topog- 
raphy, A(y) = sin(py)/3 in unit of ciq. The coupling 
strength is plotted in Fig[3]with the same parameters as 
FigH 
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FIG. 4: Absolute value of the coupling strength for the same 
2-D trilayer in Fig[3]is replotted as a function of p for a fixed 
D. The solid/gray lines are for D = 6/7 ao, respectively. The 
resonance peaks appear near p ss 1/D. 



Since the gray and dashed lines are for different p, it 
is not surprising that they oscillate with different period. 
However, when D is much larger than the amplitude of 
both topography, the roughness becomes immaterial and 
both lines start to merge with the solid line for smooth 
interfaces. This is consistent with the experimental ob- 
servation that the effect of irradiation diminishes when 
the spacer gets thicker (l2j. The important feature to 
note here is that, unlike Fig [2] where uncorrelatcd rough- 
ness always diminishes the coupling strength, the cou- 
pling here can be enhanced at certain D, as was observed 
experimentally [l^] • 

The coupling constant J is plotted in Fig[4]as a func- 
tion of p for two fixed values of D. Two signatures are 
worth noticing: (1) both lines exhibit two peaks - the 
sharp one has been discussed in section II to occur at 
p w 1/D and is caused by the bound state, while the 
broad one at p f=a kp comes from the scattering state 
where we chose kp w 1; (2) the coupling strength de- 
creases when p > 1 . 

Let us now generalize to a 3-D trilayer. Again, as- 
sume the two interfaces to share the same topography, 
A(y, z) = l/3sin(y / n y) sin(y / n z), after being irradiated 
by ion beams with flucncc n. The variation of the cou- 
pling strength with respect to n is plotted in FigfSl which 
resembles the gray line for the 2-D case in FigJU The 
position of the first peak is shifted down to such a small 
value that it is hard to discern it in the figure. Note that 
both the qualitative behavior and position of the second 
peak in Fig[5] agree excellently with the solid-circle line 
in Fig. 3 of Ref.[12| by Demokritov et al. The reason why 
we did not extend the ion dosage to higher values is that 
they violate the requirement pA <C 1 mentioned in the 
beginning of section II for the perturbation method to 
work. 
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FIG. 5: Absolute value of the coupling strength (in unit of 
10" 5 Ry I '(2tt 3 ao)) is plotted as a function of ion dosage (in 
unit of ions/eto) f° r a 3-D trilayer. Same parameters as the 
last three figures. The average spacer width is D = 7ao and 
the topography is A(y,z) = l/3sin( v / n y) sin( v / n z). The 
resonance peak appears near "n = 0.0025", consistent with 
the condition for the bound state to become resonant, and so 
is hard to discern in the figure. This peak is sharper than 
that of Figf?] for a 2-D case with a magnitude close to zero. 

IV. CONCLUSIONS AND DISCUSSIONS 

In conclusion, wc have developed a pcrturbative 
method to handle the interface roughness. The broad- 
ening of the peaks in the spectrum and the usual trend 
of suppressing the interlayer coupling can be calculated 
quantitatively by our procedures, but they are already 
captured by the previous method of static average. What 
is special about our method is that it can predict an 
enhancement in the interlayer coupling, as has been re- 
ported cxpermcntally, when the Fourier transformation 
of the topography on both interfaces is dominated by the 
same prominent component p which is comparable to kp. 
We ascribe it to the scattering states under such a cor- 
relation. In the mean time, if p is close to the inverse 
of the average spacer width, wc show that the bound 
state, namely, the peaks in the spectrum will be shifted 
and give rise to another sharper anomaly in the inter- 
layer coupling, which can be either an enhancement or 
suppression depending on the specific value of D and p. 
Under the classical mechanism of dipole interactions fl9|. 
correlations between the topography on both interfaces 
have been studied to give rise to an enhanced coupling. 
However, it is precipitated by the largeness of IR ampli- 
tude instead of when p ps 1/D or kp. 

Biquadratic coupling [il, has been known to coexist 
in the trilayers with the bilinear one that we calculated 
in this work. For completeness, we have also checked 
the effect of IR on the intrinsic biquadratic coupling by 
the perturbative approach. It is found to be enhanced 
as well under the same conditions as the bilinear term. 
However, the enhanced strength remains three orders of 
magnitude smaller than the latter, which conclusion is 
similar to that reported^ for a smooth interface and 
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without the intralayer coupling. 

Another common feature linked to the ion-beam irradi- 
ation is the alloying effect. Although this is more likely to 
take place at the interface, the possibility of its happen- 
ing in the whole spacer can not be excluded. When this 
occurs at some special composition, the nesting effect [23| 
might come about and give rise to an enhanced intcrlaycr 
coupling. The effect of alloying at the interface is consid- 
ered in Appendix D. It is found to always suppresse the 
interlayer couplin g, c onsistent with the conclusion of ab 
initio calculations [25| . When the interdiffusion is severe, 
the reduction is exponential. If the extent of alloying is 
narrow, the correction is small and the relative enhance- 
ment due to the correlations is the same as in the absence 
of the alloying. 

The excellent agreement of our conclusions and figures 
did not rely on any fitting parameter, and the micro- 
scopic topography of the interface roughness is our only 
input. We thus conclude that the enhancement of inter- 
layer coupling strength, the prime interest of this work, 
does not depend critically on the anisotropic field due to 
the interface roughness, nor the possible new channel for 
direct exchange between the magnetic side layers which 
is made possible by the irradiation [l2^. 

We thank Professors John Cheng-Chung Chi and 
Raynien Kwo and Mr. C. H. Chen for useful comments, 
and acknowledge the support by National Science Coun- 
cil in Taiwan under grants 95-2120-M007-008 and 96- 
2120-M-007-002. 



where $o(x,y) and *S>o(x,y) are the original scattering 
states for a smooth interface, and ^o(x,y) carries mo- 
mentum (k' x ,k y ). For an elastic scattering, the dispersion 
relation in Eq. (|A2[) is: 
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where m* denotes the effective mass of the charge carri- 
ers. Insert Eq. (TA"3|) into Eq. (|A"2j) . and use k x A(y), k' x A(y) 
as the perturbation to expand the boundary conditions 
in Eq. (|Al[) . We shall retain only up to the second order 
in perturbation for a ky ^ y and b ky , qy : 
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APPENDIX 

APPENDIX A: PERTURB ATIVE METHOD OF 
SCATTERING 

In order to realize the IR in a microcsopic way, we 
start from the microscopic Hamiltonian for a 2-D het- 
crojunction with an irregular topography, A(y), and the 
potential on the right side is set to be higher by Vb- The 
wavefunctions on the left/right sides are respectively de- 
noted by $>[x,y) and ^(x,y). The boundary conditions 



$(A(y),y) = V(A(y),y) 
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&tf{x,y) 



x=A(y) 



dx 



x=A(y) 
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In order to find the scattering wavefunctions, let us 
consider a plane wave with momentum (k x ,k y ) moving 
from the left towards the interface. In the limit that 
A <C l/k x , 1/p, the scattering states can be obtained by 
treating the IR as a perturbation to the smooth interface: 

*(x, y) = <t> (z, y) + J2 *k y , qy e- iq * x+i9 » v 

* (a:, y) = * (», v) + E K^ x+iqvV ( A2 ) 



where the superscripts (1/2) denotes the first/second or- 
der result, the subscript "k y ,q y " denotes scattering from 
k y to q y state, T k3 , k i is the transmission coefficient for a 
smooth interface, and u J 7 \A(y),q y — ky]" represents the 
Fourier Transformation of function A(y) to the momen- 
tum space in q y — k y . 

When !F[A(y),p\ contains many p components, their 
contributions to the third and fourth terms in Eq.{T]) 
tend to cancel each other at the level of first order per- 
turbation, unless one or two p dominants. When the 
cancellation happens, we are forced to go up to the sec- 
ond order. Except when the scattered wave happens to 
exhibit the same momentum as the incoming wave, the 
different components again cancel each other and we are 
left with the same result as that of the static average: For 
example, A(y) = J2 n [B n sm(p n ) + C n cos(p n y)] contains 
many modes, the perturbative wavefunctions in Eq. (|A2[) 
are reduced to: 
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e iKx+ik y y + R K K (\ - k 2 x d 2 )e- ik * x+ik y y 
(k x — k x ) 2 d,Q 
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where k 2 x d\ = k x J2 n (B 2 + C 2 ) signifies that they origi- 
nate from the second-order perturbation. The above re- 
sults can be shown to be equivalent to the effect of suc- 
cessive scattering from two smooth interfaces with V2do 
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of R R . The notation if ± denotes 
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V2 
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do 
V2' 
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FIG. 6: As a plane wave is scattered by a rough interface, 
the second-order effect on the wavefunction can be shown to 
be equivalent to the successive scattering from two smooth 
interfaces. 
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When the topography A^y) / A R (y) at the left /right 
interfaces is considered, the second-order perturbation 
result in Appendix A will introduce corrections to the 
reflective matrices {R~[ + + £l)/{R r ~ + sr)'- 
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and similar expression for i R after replacing the L by R. 
By use of Appendix A, the weighting b^ and d can be 
copied as: 



APPENDIX B: CORRECTION TO THE 
INTERLAYER COUPLING 

Interlayer exchange coupling in Eq.(£3]) can be written 
in a more general way by the T MatrixQ: 
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" /R R + are the reflective 



where I is the unit matrix, R~^~ 

matrices from the left /right smooth interfaces which is a 
function of energy E and are nonzero only in the diagonal 
components, and the superscript — h denotes that the 
left-going wave is turned into right-going and vice versa 
for - 



Matrix R L + takes the formQ: 



Rl + 



r L,k ± 










(B2) 

where rL,k ± / r L,q ± are the reflective coefficients of mo- 
mentums k±/qj_ from the left interface. Similar expres- 
sion by replacing the subscript L by R is the definition 



i{k±_ - k' RA _)T Rff -,T[A R {y), k y - q y ] 
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Insert the full reflective matrices with IR into Eq. ()B1[) 
and expand the trace term to the second order in sl/sr- 
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e =e L e ifc+D R+-e ifc ~ D + Rl + e ik+D e R e ik ~ D 

+ e L e ik+D e R e ik - D (B8) 

Substituting the above two definitions into Eq. (|B6[) . the 
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intcrlaycr exchange coupling can be obtained as 



take the limit Al .r(u) ~~ * or p — > 0. 
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APPENDIX C: ENERGY SPECTRUM IN ROUGH 
THIN FILMS 

The photoemission intensity 1(E) can be used to de- 
termine the energy spectrum [20] of a smooth thin film, 
which is a function of the reflective coefficients r^/ji on 
the left /right interfaces and the film width D: 



1(E) oc 
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1 -r L ,k x rR, kx e i2k ^ D +iA 
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where A is the broadening factor [20j which depends on 
the film thickness and the mean free path. Bound states 
satisfy rL,k x r R,k ± e ± = 1 and correspond to the 
peaks in 1(E). When IR are considered, more close-loops 
from successfive reflections and the action of static aver- 
aging will modify both r L / R . Eq. (|Cl[) becomes 
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where 4,js, by\ „ , and bi 



iA 



are defined in Eq. (|B5 
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It is checked that Eq. (|C2[) reduces to Eq. ljCip when we 



APPENDIX D: ALLOYING 

Besides the topography that we have discussed in de- 
tail, the rcallignmcnt of atoms at the interface due to 
ion beam irradiation should also include alloying. We es- 
timate its effect on the coupling by softening the sharp 
boundary at the interface into a smooth Eckart poten- 
tial, V(x) = Vo[l + tanh(x/a)]/2 where a depends on 
the extent of alloying. The reflective probability for this 
potential can be found in the literature [2(| as: 



cosh[27ra(fc a; — k' x )] — 1 
cosh[27ra(fc :c + k' x )\ — 1 



(Dl) 



where k x and k' x are the longitudinal momenta in cither 
side of the interface. Comparing with the sharp barrier, 
this modulation always diminishes the reflective probabil- 
ity. When a is much larger than the Fermi wavelength, 
Eq. (|Dl[) is suppressed exponentially as exp [ — 8nk x a\. 
According to Eq.([2]), the magnetic coupling strength is 
proportional to the reflective probability, and so is ex- 
pected to be reduced drastically by the alloying as con- 
cluded by the ab initio calculations [25j . 

If the interdiffusive process is weak, namely, k x a <C 
1, the alloying only introduces a small correction to the 
interlayer coupling strength: 



AE a « AEe 



- 4 (7ra) 2 k irk'p- 



(D2) 



This means that the relative enhancement due to the 
correlations is the same as before considering the alloying. 
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